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Abstract 

In this paper we apply the global isomorphism approach [V. L. Kulinskii, J. Phys. Chem. B 
114 2852 (2010)] between the Lennard- Jones fluids and Lattice Gas model to the study of the 
liquid-vapor equilibrium for the systems with the short-ranged potentials like Buckingham and the 
^ ■ Mie-potentials. The estimates for the critical point locus correlate quite well with the available 

o: 

O . numerical data. Also within the proposed approach we give the explanation for the correlation 

between the value of the second virial coefficient at the critical temperature and the particle volume 

'> 

V^ ■ found in [G. A. Vhegenthart and H. N. W. Lekkerkerker, J. Chem. Phys. 112 5364 (2000)]. 
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I. INTRODUCTION 

The Principle of the Corresponding States (PCS) is one of the driving force in searching 
the unifying description in the variety of the thermodynamic properties of the complex 
matter. From the microscopic point of view the PCS is based on simple scaling properties 
of the interaction potential which is usually assumed to be of pairwise character. Actually 
the microscopic interactions in real substances are more complex and do not conform with 
the conditions at which the PCS can be derived rigorously from the first principles |I|, Q. 
The simple scaling form of the PCS was extended to include more parameters which are 
connected with the most important properties of the interparticle interactions. For example 
K.S. Pitzer P] generalized the PCS by inclusion the acentric factor uo to account for the 
nonsphericity effects in the interaction. In connection with this a lot of model potentials 
have been invented for the description of the thermodynamic behavior of the fluid systems. 
Owing to the fast computers and sophisticated numerical algorithms now one can relatively 
easily to get the thermodynamic properties of complex matter from simulations [§. The 
detailed description of the bulk liquid-vapor coexistence as well as interface is possible for 
the model potentials in addition to usual structural measurements in real liquids. The most 
crucial property with respect to which the qualitative division of the potentials into distinct 
classes of the thermodynamic behavior has the physical meaning is the effective range of 
interaction. Indeed, it is known from the renormalization group treatment the interactions 
decaying faster than r^'y<^+^) , s > 2 — 77/^, where 77/^ ~ 0.033 is the Fisher's critical exponent 
lead to the Ising-like critical behavior, where the fluctuations play essential role. 

Recently in series of works of E. Apfelbaum and V. Vorob'ev p, |^ new correlations 
between the locus of the critical point (CP) and the line defined by the equation 

n I 

where P, T, n are the pressure, the temperature and the density correspondingly, have been 
discovered. This line is called the Zeno-line. Along this line the compressibility factor Z is 
equal to unit which means that the ideal gas equation of state P = nT is fulfilled. For the 
classical van der Waals (vsW) equation of state (EoS) [§: 

" - - a ra^ 



nb 



the relation (|l]) leads to the perfectly straight line: 



n T . 



where rib = 1/b and Tg" = a/b is the Boyle temperature in the van der Waals approxi- 
mation [Q. The relation (^ is known as the Batchinsky law [^ (see also [0). As has been 
shown by Apfelbaum and Vorob'ev this law is fulfilled with quite good accuracy for wide 
variety of substances as well as the model systems [0, 0. The authors themselves noted 
that their results can be thought in terms of some extension of the PCS because they can be 
applied to a wider group of substances in comparison with ones satisfying the corresponding 
states principle [0. Of course it is the symmetry which underlies PCS itself as well as any 
of its possible extension. The classical PCS is based on the simple scaling symmetry of the 
Hamiltonian. But the empirical facts of p, [1^] testify that the similarity of the thermody- 
namic properties has much wider applicability then one can get from the trivial rescaling of 
the parameters of interaction. 

The straightness of the Zeno-line (||) seems rather special although the deviations are 



small at least for the simple molecular liquids and for model potentials [12, 13]- As was 
shown in |jT^ (H) is equivalent to the linear density dependence of the inverse value of the 
binary correlation function at contact. It is easy to see that in low density region the line 
(0) indeed is close to the straight line. The latter intersects the T-axes at the Boyle point 
Tb determined by the condition: 

B2iTB) = 0, (3) 



where B2 is the second virial coefficient 15 



In addition to (^ the slope of the line should be determined so that it will be the tangent 
to the extrapolation of the binodal. Usually, this is fixed by the condition |13| : 



The search for the symmetry which generalizes the classical PCS and underlies the em- 



pirical findings in B O] can be done within the approach which was proposed in |16| and 



further extended in [17 . It is based on the global isomorphism between liquid- vapor part of 



the phase diagram of the molecular fiuid and the Lattice Gas (LG) or the Ising model. It is 



defined by the Hamiltonian 

H = —J N^ rii Tij — /i 2_] ^i ■ (5) 

(u> » 

Here rij = 0, 1 whether the site is empty or occupied correspondingly. The quantity J is the 
energy of the site-site interaction of the nearest sites i and j, /i is the chemical potential. 
We denote by t the temperature variable corresponding to the Hamiltonian (^. The order 
parameter is the probability of occupation of the lattice site x = (rii) and serves as the 
analog of the density. 

The aim of this work is to apply the above mentioned approach to the prediction of the 
locus of the CP for the widely used potentials and to compare the results with known . 

II. THE GENERALIZED PRINCIPLE OF CORRESPONDING STATES 



In [0 it was shown that the transformation of to the LG- variables (x, t) of the form: 

allows to map the liquid-vapor part of the phase diagram for the molecular system to the 
phase diagram of the lattice model. In particular the (|^) maps the line x = 1 which 
represents the "holeless" state of the LG onto the tangent line to the liquid branch of the 
binodal at T — )■ ||l6l . In fact the transformation (P) was proposed in [^] as the geometrical 



reformulation of these empirical facts. 

Here the parameter z can be defined by the correspondence between the loci of critical 
points of fluid {ric, Tc) and the lattice gas Xc = 1/2 , tc = 1: 

T 

(7) 



Tz-T, 

The choice of the parameters Tz and Ub will be discussed below in Section |T|. 
The coordinates of the CP for the liquid are: 

If Tz and n^ are fixed then z parameterizes the locus of the CP: 

n^ T; ^ z+l/2 

Ub Tz l + z ^ ' 



This means that for the substances belonging to the same class of the corresponding states 
represented by z, the loci of the critical points scaled to the Tz and n^, lie along the straight 
line. This correlates with the empirical analysis in 0. Note that the specific form of the 
transformation (^) and all results derived from it rely solely on the geometrical facts about 
phase diagrams of the Lattice Gas and the fluid. If the global isomorphism is accepted, (^ 
can be considered as the generalized PCS for the systems isomorphous to the Lattice Gas. 
The parameter z parameterizes the classes of corresponding states. The latter are defined 



in common way via the scaling of the parameter of the interaction |17 . 

Of course (|]) is not the exact transformation, nevertheless, as it is shown in |]l^, this 
gives the consistent framework for the description of the known empirical facts and allows 
to get rather good estimates for the locus of the CP directly from the interparticle potential. 
It is based on the simple scaling relations ric ~ ^,Tc ~ |'^att('"c)|; where r^. is the mean 
interparticle distance in the critical state. The relation 

d\r.{TjTz) 1 ^^^^ 



follows from (|]) and leads to: 

d\n 



d In ( ric/nb ) z 
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(11) 



d In ( ric/nh ) z 

If the attractive part of the potential has the power-like asymptotic $««(?") oc r^^'^^'') then 
1/z = 1 + s/d. E.g. for three dimensional LJ potential z = 1/2. The key point is the 
determination of the Zeno-line parameters n^ and Tz- Below we describe our approach to 
the determination of these parameters. 

It is still seems quite unusual that the coordinates of the the critical point which belong 



to the dense fluid region are correlated with the properties derived for low density states []10 
We believe that this fact is reflected by the global character of the proposed transformation 
(^). In the end of this section we show how the proposed approach allows to clarify the 
essence of another example of correlation between critical parameters and the low density 



properties which was discovered in [18|. It was demonstrated that the critical temperature 



for many three dimensional LJ-like systems obey the relation: 

B2{T,)^-Cvo, (12) 



with C ~ 6, fo = I cr^ and a is the distance where the potential $(r) crosses zero $(cr) = 0. 
Here 

oo 

B2{T) = 2n /" ( 1 - e-*(")/^ ) r^ dr (13) 



is the second virial coefficient. We call (|T2]) the Vliegenthart-Lekkerkerker relation. The 



results of ||T^ allow to say that the value of C for the LJ potentials appears to be close to 
2tt. Further we will use this value. 

It seems plausible that the analog of (|T^ should exist for other model systems. In such 
case the parameter C represents the class of corresponding states. Unfortunately, the exact 
connection of the value C with the specific properties of the interaction potential remains 
unclear. To clarify this issue we note that in accordance with @ for the LJ fluids the 
approximation 

Te^T^^'^^Vs. (14) 

can be used. Taking the formal high-temperature expansion of (p!3[ ) at T — ;• oo (see e.g. [§) 
we get: 

i?2(T,) = 6(l-Ti;'^VTc)+... (15) 

where all higher order terms are omitted. Since b = 4wo we get the estimate: 

B2iT,) = -8vo + ... , (16) 

which is quite satisfactory explains the value of factor C in Vliegenthart-Lekkerkerker re- 



lation (|T2|). It can be shown using the specific form of the potentials (|25D that the higher 



order terms increase the zeroth order coefficient in (|T6|). 

III. THE LOCUS OF THE CRITICAL POINT 

One of the most fascinating guesses of Apfelbaum and Vorob'ev was that the Zeno-line 
should be the tangent to the extrapolation of the binodal into the unphysical region T — )■ 
where the liquid phase ceases to exist as the stable state 0. This assumption allowed 
to consider much wider class of substances than in classical PCS [0] and formulate the 
conception of the Triangle of Liquid-Gas states H. 



It is rather tricky to extrapolate the hquid branch of the binodal into the region of sohd 
state. This is the region which determines the parameter n^. The natural parameter which 
restricts the applicability of this procedure is the locus of the triple point. 

The isomorphism transformation (^ relies heavily on the correspondence between the 
tangent to the liquid branch of the binodal and the line x = 1 which represents the "holeless" 
state of the LG |]TB[. Since any EoS can be approximated by the vdW EoS for which the 
Zeno-line is exactly the straight line we use the the generic vdW approach to the EoS 
proposed in |T^ for the short ranged potentials. It is based on the fact that the EoS for 



molecular fluids can be written as following: 

p = p+{n,T)+p^{n,T), (17) 

where p+ is the pressure contribution due to hard core repulsive interactions and p- is the 
contribution of the attractive long range part of the potential. These two contributions are 
represented in the van der Waals (vdW) form: 

Ti T 

p+ = - — , p^ = -An\ (18) 

1 — Bn 



E.g. for the potentials with the hard core the coefficients A and B are given by ||19| : 

+ 00 

A = - ^- T J drr' y{r, n, T) f'{r, n, T) , (19) 

B=b ^/"'"'^^ , , (20) 

where a is the diameter of the hard core, / = e"^"*^^'-* — 1 is the Mayer function, y(r, n, T) = 
e^'"^'^^ g{r,n,T) - the cavity function [jl5|, g{r,n,T) - pair distribution function. Thus the 



vdW-like coefficients A and B are in general the functions of the thermodynamic state 
which are in general nonanalytic |T^. Once the cavity function is determined the vdW-like 
coefficients can be obtained as the corresponding low density limiting values: 

limA= a + 0{l/T), b = Urn B , (21) 

rn-O n— >-0 

a= IT ^{r)r^dr (22) 

because if n — )■ then y{r, n,T) — )■ 1. Thus the Boyle temperature in the vdW approximation 
is: 

yHM/) ^ « _ (23) 



The parameter rib is determined by the equation analogous to the one used for the Zeno- 

(vdW) . 



hne |T0[ with the change Tb ^ Tz = T^ 






In this section we apply the mean-field scaling approach proposed in |T^ for the deter- 
mination of the CP locus. It is based on the scaling properties of the attractive part of the 
potential and the global isomorphism and gives rather good estimates for the LJ fiuids. 

We consider the set of generalized LJ potentials (so called Mie{m,n) potentials) |2^ (see 

also [|1]): 

m / ra \"/(«-™) //cr\" /o-x^x , ^. 

Hr;m,n)=e - - - - , 25 

n — m \m / \ \ r / \ r / / 

Note that at n — )■ cxd we get the Sutherland potential 

{oo r < a, 

(26) 
-;lr r> a . 

Further we pay special attention to the case ?7i = 6, for which the extensive numerical results 
are available p^-p4|. In particular, the results obtained for the liquid- vapor coexistence 
curve in p2| for the set $(r;6,n) with 7 < n < 32 testify that scaled to the critical 
parameters these curves practically coincide. Thus the PCS holds in this case. This fact 
stipulates the search of the scale symmetry which connects the critical temperatures for 
$(r;6,n) with different n. 

This can be done on the basis of the scaling property of the attractive part of the potential 



since both Tz and Tc are determined by it in accordance with |]16[. Since Tz is determined 



by the T^ there is the arbitrariness in the choice of the effective diameter of the particle. 
This scale represents the separation of the regions attributed to the repulsive and the at- 
tractive part of the potentials. Usually, it is determined by the distance where the potential 
crosses zero. For the potentials (^) it is given by a. Then the potentials $(r;6,ra) differ 
only by the scale factor at r~^ and in accordance with (p3D we obtain: 

(27) 

We choose the standard LJ potential $(r; 6, 12) as the reference case. Then from simple 
rescaling we get the dependence of the critical temperatures for $(r; 6,n): 



rr"^(6,r.i) 
Tr^^(6,n2) 


/"'2 - 

V-1- 


^ '/i-i no 


- / '^1 \ 2(ni-6) f n2 ^ 2(n2-6) 
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Uc 


Te 


^iVL) 
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0.269 


2.801 


2.107 


8 


0.303 


2.053 


1.832 


9 


0.314 


1.732 


1.648 


10 


0.319 


1.546 


1.516 


11 


0.321 


1.423 


1.417 


12 


0.322 


1.333 


1.338 


14 


0.323 


1.212 


1.223 


15 


0.322 


1.168 


1.178 


16 


0.322 


1.132 


1.141 


18 


0.32 


1.075 


1.079 


20 


0.319 


1.031 


1.032 


32 


0.308 


0.897 


0.875 



TABLE I. The results for the locus of the CP obtained from the scale invariant mean-field approach. 
The values Tc obtained using Vliegenthart-Lekkerkerker relation (^) with C = 2tt. 



It is interesting to note that both estimates (12) and (|14|) practically coincide for the standard 
LJ potential $(r; 6, 12). The relation 



T^'^^{6,n) = 3T^^^\n) 



(29) 



gives n f» 11.94, which is very close to n = 12. Note that Vliegenthart-Lekkerkerker relation 



(|T^) does not depend on the choice of the distance separating the attractive and repulsive 
regions, i.e. it is scale invariant. 



The values for the critical density ndn) are calculated according to (p^). The results are 
shown in Table fT| and Figs. |l|,^ For completeness we give the values of T^ obtained using 
the Vliegenthart-Lekkerkerker relation ([l^). 

The results for the Sutherland potential can be obtained from (pSj) by taking the limit 
n — 7- oo. Also we put Uba^ = 1 in view of natural hard core cut-off for this potential. 

The same procedure can be applied to the Buckingham potential: 



$B(r;a) 



oo 



,47^(fexp(a(l-x))-;^) 



X < ro/r„ 
X > ro/r„ 



(30) 
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FIG. 1. Critical temperature Tc of Mfe(6, n)-fluid: theory (open) and the numerical data (filled) 
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1 22, 24 1 . Also the results for Tc obtained with the help of the Vliegenthart-Lekkerkerker relation 
([T^ ) with C = 2iT are shown (crosses). 
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FIG. 2. Critical density Uc of Mie{6, n)-fluid: theory (open circles) and the numerical data (points) 
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0.286 
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TABLE II. The results for the Sutherland potential (Mze(6,n), n = oo) (|26 
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FIG. 3. Critical temperature Tc of Buckingham potential: theory (dashed) and the numerical data 



(crosses) pq], the Vliegenthart-Lekkerkerker relation (12) with C = 27r (stars). 



where r^ determines the hardcore cut-off and is defined by the point of maximum of $b. 
The distance r^ is the point where ^b reaches minimum $s(ro) = — £. Also the distance 
where $b crosses zero we denote a. It depends on the parameter a and clearly tq < cr < r^. 
Comparing (pSJ) for m=6 and (0) we see that the corresponding change of scale is needed. It 
is easy to find that at a ~ 12.8 the attractive part of the Buckingham potential equals to that 
of the standard LJ 6-12 potential $(r; 6, 12) and therefore for such potential T^ = T^ /3. 
So we choose this case as the reference point. The analog of (^8]) is: 



TM) = -( — 
a 



v/l-6/a 



(31) 



The density made dimensionless in units of the distance where the potential has minimum 
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FIG. 4. Critical density ric of Buckingham potential: theory (filled) and the numerical data (open) 
[25|. The dashed line is the interpolation of the calculated values of ric- 



The comparison of the results and available numerical data for the Buckingham potential 
is given in Figs. |^,|. 

IV. CONCLUSIONS 



We considered the application of the scale invariant mean-field approach as the conse- 
quence of the global isomorphism between LG and LJ-like fluids to the calculation of the CP 
locus for such fluids. Simple scaling expressions for the critical temperatures are obtained 
for Mie{6,n) and Buckingham potentials. The comparison with the available numerical 



data shows satisfactory agreement. Also we proposed the explanation for the result of [[18 



namely, the constancy of the value of the second virial coefficient at the critical temperature. 



The presented results are based on the 



asymptotical behavior of the attractive part of 



the potential. For completeness it would be interesting to apply the proposed approach to 
the non power-like potentials like Yukawa or step-like potentials, e.g. square-well potential. 
We leave this for subsequent investigation. 

In previous consideration we assumed that the liquid-vapor binodal as well as the CP 
exists. But it is known both from the computer simulations and theoretical models that 



the region of the liquid phase may collapse if the range of interaction is too short (see |T^ 



12 



and reference therein). This happen if the triple point meets the critical one. In such 
situation there is no critical state because the liquid branch is not the equilibrium phase in 
this situation. 

Obviously the sufficient condition for the liquid branch to disappear as a stable state is 
Tc < Tfr, where Ttr is the triple point temperature. As the value of z decreases with the 
decreasing the range of interaction the situation when Tc occurs below Ttr becomes possible. 
Therefore for the potential with varying range of interaction, e.g. via the power of the 
attractive part, there is the value z* below which the liquid branch of the binodal ceases to 
exist as the stable phase. 

From the point of view of the global isomorphism with the LG this could be interpreted 
as the moment when the interaction range becomes smaller than the lattice spacing of the 
effective LG model. 

To estimate the value of z* we use the standard result that for LJ fluid Ttr ~ 0.75 Tc ~ 
T^ /A. Thus from Ttr = Tc we get the estimate z* ~ 1/4. Thus from ( [TT| ) we get: 

rfln($(n-i/")/rz) 
d Innc/rih 

Applying the condition ( ^2|) to the Sutherland potential (|26|) and the Yukawa potential: 

oo , if r < cr 

Hr) = { (33) 

— £ ^exp ( — A(r/o- — 1) ) , if r>a, 

we get that the liquid branch loses its stability if s > 9 and A > 11 correspondingly. The 



computer simulations give the values s > 5 and A > 6 [^, ^ . Obviously, the disagreement 
is caused by the crudeness of the estimate for the triple point temperature which was taken 
as for standard L J fluid. In fact the position of the triple point also depends on the potential. 
The computer data correspond to the estimate z* ~ 1/3. This issue will be the subject of 
the separate paper. 
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